NASA/TM— 2000-210057 



A Wideband Satcom Based Avionics Network 
With CDMA Uplink and TDM Downlink 

D. Agrawal and B.S. Johnson 
University of Illinois, Urbana, Illinois 

U. Madhow 

University of California, Santa Barbara, Santa Barbara, California 
K. Ramchandran 

University of California, Berkeley, Berkeley, California 
K.S. Chun 

Glenn Research Center, Cleveland, Ohio 


4 


June 2000 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 

has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or co sponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized sendees that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
data bases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http:llxvivw.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at (301) 621-0134 

• Telephone the NASA Access Help Desk at 
(301) 621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 


NASA/TM— 2000-210057 



A Wideband Satcom Based Avionics Network 
With CDMA Uplink and TDM Downlink 

D. Agrawal and B.S. Johnson 
University of Illinois, Urbana, Illinois 

U. Madhow 

University of California, Santa Barbara, Santa Barbara, California 
K. Ramchandran 

University of California, Berkeley, Berkeley, California 
K.S. Chun 

Glenn Research Center, Cleveland, Ohio 


Prepared for the 

18th International Communications Satellite Systems Conference and Exhibit 
sponsored by the American Institute of Aeronautics and Astronautics 
Oakland, California, April 10-14, 2000 


National Aeronautics and 
Space Administration 


Glenn Research Center 


June 2000 



This report contains preliminary 
findings, subject to revision as 
analysis proceeds. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 
Price Code: A03 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 
Price Code: A03 




A Wideband Satcom Based Avionics Network With CDMA 

Uplink and TDM Downlink 


D. Agrawal and B.S. Johnson 
University of Illinois 
Urbana, Illinois 

K. Ramchandran 
University of California, Berkeley 
Berkeley, California 


Abstract 

The purpose of this paper is to describe some key 
technical ideas behind our vision of a future satcom 
based digital communication network for avionics 
applications The key features of our design are as 
follows: 

(a) Packetized transmission to permit efficient use 
of system resources for multimedia traffic; 

(b) A time division multiplexed (TDM) satellite 
downlink whose physical layer is designed to oper- 
ate the satellite link at maximum power efficiency. 
We show how powerful turbo codes (invented orig- 
inally for linear modulation) can be used with non- 
linear constant envelope modulation, thus permit- 
ting the satellite amplifier to operate in a power 
efficient nonlinear regime. 

(b) A code division multiple access (CDMA) satel- 
lite uplink, which permits efficient access to the 
satellite from multiple asynchronous users. Closed 
loop power control is difficult for bursty packetized 
traffic, especially given the large round trip delay 
to the satellite. We show how adaptive interference 
suppression techniques can be used to deal with the 
ensuing near-far problem. 

(c) Joint source-channel coding techniques are re- 
quired both at the physical and the data transport 
layer to optimize the end-to-end performance. We 
describe a novel approach to multiple description 
image encoding at the data transport layer in this 
paper. 

1 Introduction 

An upgrade of the telecommunications infrastruc- 
ture for avionics applications will require the use 
of advanced digital communication methods which 
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support a diverse mix of traffic requirements, in- 
cluding voice, data, image, and video transmission. 
Such a system w r ould involve not only communica- 
tion between aircraft and ground stations, but also 
between aircraft and aircraft. Satellite communica- 
tion forms an attractive option for providing global 
connectivity in this scenario. This paper reports 
on a research effort aimed at designing the physical 
layer for a satcom based digital avionics network. 
We propose to utilize the Ku and Ka bands. Since 
the link to and from the satellite is expected to form 
the bottleneck, we focus on the design of that link 
in this paper. In order to minimize the interference 
caused to other satcom systems in the same band, 
spread spectrum transmission will be employed, so 
that the transmitted power is spread over the entire 
available bandwidth. Furthermore, in order to fa- 
cilitate support of multimedia traffic, packet-based 
transmission will be used. Some of the salient fea- 
tures of the design, and of the results presented in 
the paper, are summarized below'. 

Downlink: The one-to-many “dowmlink” from 
the satellite is severely powerlimited, with the result 
that the power amplifiers in the satellite typically 
operate in a nonlinear saturation mode. Hence, 
constant envelope modulation will be employed on 
the downlink. Continuous phase modulation (CPM), 
w-hich is a nonlinear modulation strategy, is a band- 
width efficient strategy in this context. However, 
it has not become as popular as linear modulation 
strategies, because the state of the art of coding 
techniques for CPM is not advanced. In Section 
2 of this paper, w r e show r how' powerful turbo-like 
codes, developed originally for binary linear modu- 
lation, can be leveraged into CPM systems. 

Uplink: Direct sequence (DS) code division mul- 
tiple access (CDMA) with linear data modulation 
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will be employed for the many-to-one “uplink” to 
the satellite. The data for each user modulates 
a user-specific wideband spreading w'aveform. We 
employ short spreading waveforms, in which the 
period of the spreading waveform equals the sym- 
bol interval. This is in contrast to the CDMA- 
based US digital cellular standard, IS-95, which 
employs “long” spreading waveforms whose periods 
are much longer than the symbol interval. Short 
spreading sequences are chosen because they are 
compatible with adaptive interference suppression 
schemes [1, 2, 3, 4] that provide large gains in per- 
formance and capacity compared to conventional 
systems based on long spreading waveforms. In Sec- 
tion 3 of this paper, we provide a tutorial on recent 
developments in adaptive interference suppression, 
and give new results on joint Doppler estimation 
and interference suppresion. 

Joint Source- Channel Coding : Separate source 
coding and channel coding is optimal in an information- 
theoretic sense (i.c., with no delay or complexity 
constraints) for singleuser transmission over a time- 
invariant channel. However, joint source-channel 
coding is important for optimizing end-to-end Qual- 
ity of Service in many practical settings, including 
real-time traffic, multiple users, and time-varying 
channels. This is especially crucial for wireless im- 
age and video transmission. In Section 4 of this pa- 
per, we describe a source-channel coding strategy 
for end-to-end quality of service on a packetized 
network. Ongoing work on source-channel coding 
for a single satellite link is not described. 

After describing the main technical results in 
Sections 2 through 4, we give our conclusions in 
Section 5. 


2 Forward Link: Turbo Coded 
CPM 

Continuous phase modulation (CPM) signaling pro- 
vides bandwidth-efficient constant envelope mod- 
ulation schemes that suitable for downlinks from 
satellites. The excess phase of such signals varies 
continuously and is given by the convolution of an 
Af-ary data sequence {a*} with a phase response 
function q(t). Thus, a CPM signal can be described 


by 

s(i) = y^cos(27r/o< + 2irh £" =0 
nT ^ t < (n + 1)T 

( 1 ) 


where h = K/P is the modulation index and T 
is the symbol duration. Usually, M — 2 k for an 
integer k , and ajs take values in the set {±(Af — 
1 ),±(M - 3). ... , =L1}. Denote the excess phase of 
the signal given above by 4>(t ), that is, 

71 

< p(t ) = 27t/* ^ c*iq(t — iT) nT ^ t < (n • + 1)T 
?— o 

At present, some of the best coding schemes for con- 
tinuous phase modulation use convolutional codes 
based on Z p, the ring of integers modulo P. For 
good performance these convolutional codes have to 
be judiciously chosen for each specific CPM signal- 
ing according to the parameters h and A/, and the 
phase response function q(t) [5]. Furthermore, the 
choice of codes is greatly influenced by the channel 
and receiver characteristics. 

Recently discovered graph-based codes (most no- 
tably turbo codes) and associated iterative decod- 
ing algorithms have demonstrated coding gains that 
far exceed the coding gains provided by convolu- 
tional codes. In fact, for a wide class of binary-input 
channels and receivers, these codes can achieve per- 
formance close to the capacity [6]. However, graph- 
based codes are usually compatible with binary- 
input channels, and they cannot directly be ap- 
plied to an Af-ary input continuous phase modu- 
lated (CPM) channel. Arguably, it is possible to 
extend binary graph-based codes to Af-ary alpha- 
bets in order to use them with CPM. This approach, 
however, will not be able to capitalize on the vast 
research being done on binary graph-based codes. 
More importantly, such an approach may be un- 
necessary as long as an Af = 2* -ary input CPM 
channel can be divided into parallel binary-input 
channels. 

In this paper, we use bit interleaved coded mod- 
ulation (BICM) to convert an Af-ary CPM channel 
into several binary-input channels [7], and thus en- 
able graph-based coded modulation on CPM chan- 
nels. In Section 2.1, we briefly describe a BICM 
system model that uses CPM channels. Next, in 
Section 2.2, we describe the results of our extensive 
simulations. 
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2.1 System Model 

Figure 1 shows the block diagram of a communica- 
tion system which uses bit interleaved coded mod- 
ulation with a CPM channel. In this setup, the 
binary codeword at the output of the channel en- 
coder is fed to a bit interleaver 7r. The output of 
bit interleaver tt is then passed to a (one-to-one) 
labeling mapper /i which maps a group of A r k bits 
to a sequence of N symbol in an M = 2* -ary in- 
put alphabet X . Denote the resulting sequence of 
length A r by x. The sequence x is transmitted using 
continuous phase modulation. Note that due to the 
intersymbol interference introduced by the phase 
response function g(f), a continuous phase modula- 
tor is not memoryless, rather, it can be modeled as 
a finite state machine. We assume that in the be- 
ginning of transmission (of each sequence) the mod- 
ulator state is zero, and at the end of each sequence, 
some extra symbols are sent to reset the modulator 
state back to zero. 

At the receiver, the received signal y is passed 
to a demodulator which computes the bitwise soft 
information p(y]/ l (x) = b). Here P(x) denotes the 
i - th bit of /i _1 (x) and b belongs to the set {0,1}. 
Finally, the bitwise soft information is deinterleaved 
and is passed to a turbo decoder. Using the anal- 
ysis of BICM by Caire, Taricco, and Biglieri [7], it 
is clear that, under the assumption of ideal inter- 
leaving, the above arrangement converts an M -ary 
input CPM channel into Nk parallel, independent, 
memory less, binary-input channels. These channels 
have roughly the same SNR and noise characteris- 
tics. 

Note that, apart from the CPM demodulator, 
this setup uses standard off-the-shelf components, 
and thus reduces the implementation complexity. 
In the following subsection, we will briefly describe 
an easily implementable CPM demodulator which 
can generate bitwise soft information. This demod- 
ulator is a straight-forward result of combining Ri- 
moldi’s time- in variant trellis of CPM signals with 
theBahl-Coke-Jelinek-Raviv(BCJR) algorithm. For 
more details on these, the reader is referred to [8], 
[9], 

2.1.1 CPM demodulator with bitwise soft 
decisions 

It is well knowm that continuous phase modulated 
signals given by (1) can be represented by a time- 
variant trellis that describes the evolution of the 


{Turbo Encodcrf 
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Modulator 



Demodulator 


Figure 1: A turbo-coded and continuous phase 

modulated communication system 


excess phase <fr(t). In principle, we can obtain the 
bitwdse soft information by applying the BCJR al- 
gorithm on this trellis with a complexity that grows 
only linearly in the number of bits [9] . However, the 
time- variant nature of the trellis makes the termi- 
nation of the trellis in the zero state at the end 
of Nk bits difficult and the large number of states 
2PM l_1 in such a trellis causes extra computation 
complexity. Fortunately, an alternative representa- 
tion of continuous phase modulation, proposed by 
Rimoldi [8], results in a time-invariant trellis with 
lesser number of states. 

In Rimoldi’s formulation, instead of using (1) 
to construct phase trellis, an alternative represen- 
tation of the CPM signal is used [8]. In this rep- 
resentation, the center frequency is chosen to be 
f i — f 0 — h(M — 1) /(2 T) so that the excess phase of 
a CPM signal becomes, for t = nT + r (0 ^ r < T), 

rp(t = nT + t) = 27r/i <**<?(* ~ *T) 4* 

= 2*h Er=o L Ui + EfrTo 1 U n -if(T + iT) + W(t) 

( 2 ) 

where 

W(t) = _ 2nh(M - 1) Ef^ 1 f(r + iT) 

+ (L- 1 )(M - l)nh 

represents data independent terms, and Ui = (a l + 
(M —1l))/2 take integer values between 0 and M — 1. 
It was shown in [8] that the excess phase given 
by (2) can be represented on a time-invariant trel- 
lis with PM L ~ l states. Since this trellis is time- 
invariant and includes 0 in its symbol alphabet, it 
can be easily terminated in the zero state. 

At the receiver, the received signal is lowpass fil- 
tered and is sampled every T /£ second, where £ is 
the oversampling factor. These samples are passed 
to a CPM demodulator which computes soft infor- 
mation of the transmitted bits by using the stan- 
dard BCJR algorithm on the time-invariant trellis 
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described above. For more details on branch metric 
calculations and the effect of parameter £, we refer 
the reader to [10]. 

2.2 Simulation Issues and Results 

In order to demonstrate the performance of CPM 
modulated channels, we used classical turbo codes 
with identical constituent convolutional codes gen- 
erated by the polynomials {D 4 + D 3 + D 2 + D } -f 
1 7 D 4 + 1}. The parity bits were punctured alter- 
natively to generate a rate 1/2 turbo code. The 
interleaver length was 2 15 = 32 768. Turbo code 
with the aforementioned parameters has been eval- 
uated on a variety of channels and it serves at a 
good performance benchmark. 

A CPM signaling scheme can be specified by 
three parameters: number of data symbols M , the 
frequency response pulse g(t) = dq(t)/dt and the 
modulation index ft. We tested the performance 
of CPM channels for M = 2.4 and 8. In this pa- 
per, we provide results only for M — 4 and 8. It 
is well known that for M = 2, CPM is an inferior 
modulation scheme when compared to the CPM for 
M — 4 and 8 [10, page 182]. Our simulation results 
affirm this fact for turbo-coded CPM channels also. 

For full response signaling, we simulated contin- 
uous phase frequency shift keying (CPFSK) which 
is a special case of CPM modulation with rectangu- 
lar frequency pulses. For partial response signaling, 
we used raised cosine frequency pulses of duration 
3 T (also referred to as 3RC CPM signaling), where 
T is the symbol period [10, page 52]. 

The choice of modulation index ft = KjP is 
more involved. For a given M and g(t ), increas- 
ing ft results in larger occupied bandwidth. How- 
ever, increasing ft up to a certain limit, usually also 
results in larger minimum distance between CPM 
signals. Since the minimum distance dictates un- 
coded bit-error rate at high SNR’s, by increasing ft , 
one can trade-off bandwidth with transmitted en- 
ergy (at least for high SNR’s). After a certain limit, 
there is no benefit in increasing ft since that does 
not yield larger minimum distance [10, Chapter 3]. 

However, the trend of increasing minimum dis- 
tance with ft does not hold for certain values of ft, 
called w r eak modulation indices. . Weak modulation 
indices generate CPM signals whose excess phases 
differ by a multiple of 27T and therefore are identi- 
cal to each other. This results in low minimum dis- 
tance. Such weak modulation indices are usually 


M 

2 

4 

8 

9(t) 




CPFSK 

7 

10 

9 

10 

9 

10 

3RC 

4 

5 

1 

1 


Table 1: Optimum values of ft for different CPM 
parameters 


avoided in practice since the corresponding CPM 
schemes have low asymptotic gains. 

For turbo-coded CPM signaling, the trade-offs 
described above are invalid because turbo-codes work 
at low SNR’s, where apart from the minimum dis- 
tance between CPM signals, second order effects 
are also as important. As a rule of thumb, for 
turbo codes to work well, the uncoded bit error rate 
(BER) should be between 10-15%. Thus, w T e are 
interested in the SNR required by CPM schemes 
(with different ft) to achieve a bit error rate of 13%. 
Moreover, we are also interested in the values of ft 
that have low P } since the complexity of demodula- 
tor is directly proportional to P. In order to choose 
ft, therefore, we simulated uncoded CPM for dif- 
ferent values of ft, and chose an ft that achieved 
a good trade-off between the value of P and the 
SNR required to achieve a 13% bit error rate . For 
M = 2,4 and 8, and for CPFSK and 3RC, these 
values are given in Table 2.2. Note that for 3RC 
scheme, we have chosen weak modulation indices 
for M = 4 and 8. 

Figures 2-5 show the simulated performance of 
turbo-coded CPM signaling. These figures also show 7 
the corresponding performances of turbo-coded M - 
ary PSK modulation. Note that, w T hile PSK mod- 
ulation produces a constant envelope signal in dis- 
crete time, the analog signal corresponding to PSK 
modulation, when used with a bandwidth efficient 
pulse such as a raised cosine, is not constant en- 
velope. Thus, the satellite amplifier would need to 
back off from the saturation regime to support PSK, 
leading to a power penalty that is not accounted for 
here. Figure 2 shows that for M = 4, CPM takes 
about 1.6 dB more SNR to achieve a bit-error rate 
of 10“ 5 . However, this would be more than compen- 
sated by the ability to operate the satellite amplifier 
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Figure 2: Turbo-coded CPFSK with M = 4 versus 
turbo-coded QPSK 

in the nonlinear regime for CPM. Furthermore, note 
that power spectrum of CPM scheme is much more 
compact than the QPSK scheme. The power spec- 
trum of the CPM signaling can be made even more 
compact without increasing transmitted power by 
using the 3RC CPM signaling (see Figure 3). 

Even without accounting for the satellite ampli- 
fier operating regime, the penalty in transmitted 
power with respect to PSK modulation virtually 
disappears for M = 8. Figures 4 and 5 show that in 
order to achieve a bit-error rate of 10“ 5 , both CPM 
and PSK signaling use about the same transmitted 
power. However, the use of CPM signaling results 
in a much better power spectrum characteristics. 
We note that for M = 8, the number of states in 
the CPM trellis is 10 for CPFSK and is 64 for 3RC. 
On the other hand, the spectral characteristics of 
3RC CPM signaling are much better than those of 
CPFSK signaling [10, page 147] Thus 3RC CPM 
signaling provides better spectrum at the cost of 
more complexity in the receiver. 

3 Reverse Link: Asynchronous 
CDMA with Interference 
Suppression 

In a DS-CDMA system, the information-bearing 
signal for each user is spread over a wider band- 
width by multiplying by a spreading waveform unique 
to that user. Multiple access capability is obtained 



Slgnal-to-nolso ratio (dB) 


Figure 3: Turbo-coded 3RC CPM with M = 4 ver- 
sus turbo-coded QPSK 



Slgnal-to-nois© ratio (dB) 

Figure 4: Turbo-coded CPFSK with M = 8 versus 
turbo-coded 8-PSK 
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Figure 5: Turbo-coded 3RC CPM with M = 8 ver- 
sus turbo-coded 8-PSK 

by designing the spreading waveforms to have small 
crosscorrelations. Conventional DS-CDMA systems 
(as in IS-95) use spreading waveforms whose period 
is much larger than the symbol duration , and rely 
on interference averaging to provide robust perfor- 
mance. However, if the spreading waveforms have 
period equal to the symbol interval (these are termed 
short spreading codes), the multiple-access inter- 
ference (MAI) becomes cyclostationary, and it is 
possible in theory to obtain much better perfor- 
mance than for conventional systems by learning 
the structure of interference, i.e., by using adaptive 
interference suppression. However, a key bottle- 
neck to realizing the potential gains from adaptive 
interference suppression is that conventional adap- 
tive algorithms cannot cope with the fast frequency- 
selective fading typical of wireless channels. How- 
ever, we have recently obtained algorithms, based 
on an invention termed the differential Minimum 
Mean Squared Error (DMMSE) criterion [4], that 
are robust over wireless channels. With this recent 
invention, a practical system architecture based on 
adaptive interference suppression is within reach. 
This technology also combines naturally with an- 
tenna arrays (using either beamforming or diver- 
sity techniques) to further enhance capacity and 
performance. Finally, the use of linear interference 
suppression reduces the multiuser CDMA channel 
approximately to a singleuser AWGN channel with 
fading, which means that work on the hot topic of 
source and channel coding for singleuser wireless 
channels can be leveraged into our application. 


The new research reported in this paper con- 
sists of extending the DMMSE criterion to pro- 
vide carrier frequency acquisition as w^ell as inter- 
ference suppression using digital signal processing, 
after bringing the received signal dowrn to baseband 
using a fixed analog carrier generator. The satel- 
lite uplink is well modeled as a multiuser Additive 
White Gaussian Noise (AWGN) channel, with each 
user undergoing a (possibly different) Doppler shift. 
However, we present algorithms that can cope not 
only w T ith a Doppler shift, but also with variations 
in amplitude and phase due to fading. This is so 
that the set of algorithms developed are applica- 
ble to terrestrial wireless settings as well, as might 
be the case in communication between aircraft and 
control towers. 

Frequency acquisition consists of estimating and 
correcting for the Doppler offset suffered by the user 
of interest. We employ the DMMSE algorithm to 
provide interference suppression together with fre- 
quency acquisition. If the channel does not exhibit 
significant fading, then DMMSE adaptation can be 
replaced by standard MMSE adaptation. 

In this section, we describe algorithms based on 
the DMMSE criterion for a general system model, 
of which asynchronous CDMA with short spread- 
ing sequences is a special case. This is followed 
by a description of the current work on extending 
these ideas to provide frequency acquisition as w r ell 
as interference suppression. 

3,1 System Model for Adaptive In- 
terference Suppression 

It is useful to consider a generic model for inter- 
ference suppression (which applies to other con- 
texts such as equalization or beamforming as well) 
described in complex baseband, discrete time, no- 
tation. We will later describe how asynchronous 
CDMA with short spreading sequences is a special 
case of this model. 

It is assumed that the received signal has al- 
ready been passed through an appropriate filter and 
sampler. A linear receiver c operates on a sequence 
of received vectors {r n }. Each received vector r n is 
an L - dimensional complex vector w^hich is a sum of 
a number of signal vectors, plus noise, w'here each 
signal vector corresponds to ISI or MAI. The out- 
put at time n of the linear receiver c to be designed 
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is given by < c,r n >= Ya=i c *M r »[* ], where 

K 

r„ = B x [n]ui + B k [n]u k + w„. (3) 

k = 2 

Here ui is the desired signal vector, and u l5 ...,uj 
are interference vectors. These are assumed to re- 
main roughly the same (or to vary slowly) with n. 
The terms {B*[n],fc = 1,..,A'} can vary arbitrar- 
ily with n, and their interpretation depends on the 
application (see Sections 4.1.3 and 4.1.4 for two spe- 
cific applications). The vector w n is discrete-time 
noise. 

The standard linear MMSE receiver minimizes 
the Mean Squared Error (MSE) between the re- 
ceiver output and the desired sequence B\[n] : given 
by E[(< c,r n > -Bi[n}) 2 ]. Such a receiver can 
be implemented adaptively using a number of well 
known algorithms [11]. The adaptive algorithms 
are initialized by a known training sequence {t[n]} 
. Assuming that the desired sequence {£h[n]} satis- 
fies Bi[n] = g[n]t[n], where g[n] is an unknown gain, 
and that g[n] varies slowly with n, standard adap- 
tive algorithms based on the cost function E[(< 
c,r n > — t[n]) 2 } can automatically compensate for 
the unknown gain g[n]. However, this procedure be- 
comes difficult or impossible to implement when the 
variations in the unknown gain g[n] are too rapid 
for the adaptive algorithm to compensate for them. 

3.2 Differential MMSE Detection 

In the applications of interest, while the desired se- 
quence {l?i[n]} includes an unknown gain g[n], the 
gain g[n] is roughly constant over two observation 
intervals. Thus, we assume that A\[n] = 
is either known (because of training) or can be es- 
timated. The key aspect of our scheme is that we 
choose the correlator c to satisfy the following dif- 
ferential MMSE criterion: 

Minimize 

£[(< c,r„ > —A\ [n] < c,r„_i >) 2 ] (4) 

subject to 

E[|<c,r n >| 2 ] = c"Rc = l (5) 

where R — E[ r„r^], and where x /f denotes the 
conjugate transposed of x. 

It can be shown that, under standard assump- 
tions, the solution to the above problem is the same 


as that for the conventional MMSE problem. How- 
ever, our reformulation enables us to avoid tracking 
the unknown gain in the B\ [n] by using B± [n — 1] as 
a reference. Standard optimization techniques can 
be used to obtain a number of adaptive algorithms, 
a few of which are stated below. 

3.3 DMMSE-Based Algorithms 

The solution to the optimization problem (4)-(5) is 
the generalized eigenvector with the largest eigen- 
value for the following generalized eigenvalue prob- 
lem: 

Be = i/Rc 

where B = E[A* l [n]r„r%_ 1 + .4, [n]r„_irjf]. 

Block Algorithms 

Replacement of statistical expectations by em- 
pirical averages over M observation intervals leads 
to block algorithms. This is analogous to the block 
LS algorithm for standard MMSE adaptation. 

1. Compute R = £nli r n r " and B = £«=i A i M 
r„r^L x + jdifnjrn-ir^ . 

2. Use any standard algorithm to find the largest 
generalized eigenvector c for the pair (B,R). One 
possible algorithm is the power algorithm [12] ap- 
plied to R _1 B. 

Recursive Algorithms 

Instead of averaging over a block of data, the 
averages can be performed recursively, using an ex- 
ponential forget factor 0 < A ^ 1 . This is analogous 
to the RLS algorithm for standard MMSE adapta- 
tion. 

1. Update R„ = AR„_j + (1 - A)r„r^ and B n = 
AB^-i + (1 - A)[.4j [n]r„r^_j + 4i[n]r n _ir*. 

2. Update the largest generalized eigenvector c n 
for the pair (B n ,R n ). One possibility is to apply 
one or more steps of the power algorithm [12] on 
the matrix R“ l B n , where R” 1 is updated using a 
standard RLS mechanism [11]. 

Alternatives to the preceding are obtained by 
updating c using a stochastic gradient for the cost 
function (4), together with a normalization that 
projects c approximately onto the constraint sur- 
face given by (5). This is analogous to the LMS or 
normalized LMS (NLMS) algorithm. 

Cn— i 4“ < c n _! , r n _j j4 x [n] r n > 

c n 7 ^ : 

y c ^-l^- c n — 1 

( 6 ) 
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where 


A/ n = AM„_i+(l-A)r"r n (7) 

estimates the power of the received vector, and and 
where R is an estimate of the correlation matrix R. 

3.4 Application to CDMA with short 
sequences 

For a CDMA system in which user k uses spread- 
ing waveform Sk{t) to modulate all its transmitted 
symbols, the transmitted baseband signal for user k 
is given by bk y n$k{t-nT), where b kyU is the nth 
symbol of the kth user. This model applies both 
to narrowband systems with ISI and co-channel in- 
terference, and to wideband CDMA systems with 
short spreading sequences. For simplicity of ex- 
position, we restrict attention to frequency non- 
selective fading in this report. Assuming that the 
fading is slow relative to the symbol interval, the 
received signal for the kth user may be modeled as 
* 9k)(t - nT), where {F fc [n]} is a 
sequence of complex numbers modeling amplitude 
and phase fluctuations due to fading, and g k (t) is a 
normalized channel impulse response which is well 
modeled as time invariant (on the time scale of the 
receiver adaptation) . 

After filtering (typically using a chip matched 
filter) and sampling (typically at the chip rate or 
a multiple thereof), the received signal is of the 
form given by (3), where the time index n corre- 
sponds to the symbol time, the desired sequence 
Bi[n] ~ Fi[n]bi tn7 and with the desired vector ui 
denoting the sampled and windowed (by the fil- 
ter length) version of the cascade of the spreading 
waveform and the channel seen by the desired user. 
For k = 2 , ...,/C, the sequences {B k [n}} and the 
vectors u* represent contributions due to ISI and 
MAI. See [3] for a tutorial description of the mod- 
eling involved, and see [1], [13], [2] for other papers 
using this equivalent synchronous model. 

Assuming that the fading is roughly constant 
over two successive symbols (i.e., Fy [n] « Fi[n— 1]), 
we have that B i [n] « Ay [n] By [n — 1] , where Ay [n] = 
bi, n /bi r n-i is the sequence of symbols used to train 
the DMMSE detector, the desired user. Algo- 
rithms based on the differential MMSE criterion can 
now be directly applied, using an initial training 
sequence for {by in } 7 followed by decision-directed 
adaptation. The DMMSE detector leaves a phase 


ambiguity at the output, that can be resolved ei- 
ther by differential demodulation or by the use of 
pilot symbols. In standard differential demodula- 
tion, decisions on Ay\n] = ay >n would typically be 
based on the metric < c.r n >< c,r n _i >*. How- 
ever, it is also possible to get better performance for 
channels with coherence times of several symbols 
by using block noncoherent demodulation. A the- 
ory of signal design for such channels has recently 
been developed by Prof. Madhow [14], and will be 
leveraged into the current project as appropriate. 

The receiver structure is depicted in Figure 6. 


(Differential 
MMSE algorithm) 


Receive Filter — H Sampler] 


Adaptive 

Receiver 


(at multiple 
of symbol rate) 


■*- { Sampler] — 
(at symbol 
rate) 


i Differential 
Demodulator 


1 ,n a 1 ,n (decision-direc ted mode) 
(training mode) 


Figure 6: Receiver structure for interference sup- 
pression for CDMA with short spreading sequences. 
The receiver operates on samples at a multiple of 
the symbol rate. 


3.5 Frequency acquisition 

A frequency offset of Sf is manifested as a factor 
of in complex baseband. Over a symbol in- 

terval of length T s , this offset causes a phase shift 
e ?2n6fT a _ e j 5 B f ram ework of the general 

model (3), this means that 

By [n] = by ^ nS6 


Thus, the ratio 


Ai[n] ~ 


B i [re] 
By [n — 1] 


_ J60 

by^n-l 


needed for DMMSE adaptation training the cor- 
relator c requires knowledge of 56 (which in turn 
requires knowledge of the frequency offset 5f). 

For an unknown 50, our strategy is to hypothe- 
size several different values of 50 , and to run DMMSE 
adaptations in parallel for each hypothesized value. 
The best value will then be chosen based on re- 
ceiver statistics. In the context of the block algo- 
rithms described in Section 3.3, this corresponds to 
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computing a different matrix B for each possible 
hypothesized 50. However, computation of the cor- 
relation matrix R and its inverse, which is often the 
most expensive step, is the same as that required for 
a situation in which the frequency offset is known. 
Moreover, this correlation matrix computation is 
the same for all users being received, while the B 
computation is user-specific. 

We use computer simulations to verify the per- 
formance of the preceding hypothesis based frequency 
acquisition scheme proposed in previous chapters. 
We consider, for simplicity of illustration, a syn- 
chronous CDMA system with five active users, each 
transmitting a binary phase-shift keying (BPSK) 
modulated signal. The processing gain (he., the 
number of chips per bit) is 10, and short spreading 
sequences are employed. The value of the carrier 
frequency offset for each user is one-fifth of the data 
rate. 

We consider a packet ized system with a packet 
containing 500 bits: 50 training bits and 450 data 
bits. No decision-directed updates are performed in 
the simulations. Rather, training occurs at the be- 
ginning of each packet and values “learned” during 
training are used for demodulation throughout the 
duration of the packet. 

Simulations of the (uncoded) BER are performed 
with a signal-to-noise ratio (^) for the desired user 
of 20 dB. Since closed loop power control on a satel- 
lite link with a large round trip delay may be dif- 
ficult, we evaluate our design under a severe near- 
far problem, with tw T o inter ferers transmitting at a 
power 20 dB stronger than the desired user, two 
interferers transmitting at the same power as the 
desired user. In Figure 7, we plot, for each value of 
the BER averaged over 1000 packets. It is clear 
from the figure that the performance of the receiver 
is insensitive to the near-far problem, and that the 
10% uncoded BER required for good turbo code 
performance can be achieved at ^ smaller than 
5 dB. 

4 Joint Source-Channel 
Coding 

In this paper, our focus is on source- channel coding 
for image data. Instead of considering the physical 
layer, as done in the previous two sections, we con- 
sider source-channel coding at the data transport 
layer, for end-to-end quality over a packetized net- 


work. Due to space restrictions, numerical results 
are omitted, and only a sketch of the basic ideas is 
provided. See [15] for details. Joint source-channel 
coding strategies for the physical layer of a satellite 
link is a subject of our ongoing research, and will 
be reported on in future publications. 

Typical image/video bitstreams are highly struc- 
tured i.e., characterized by a natural hierarchy of 
importance layers or resolutions. In fact, there is 
a strong push in the evolving next-generation im- 
age compression standard, JPEG- 2000, to support 
bit stream scalability (bit stream scalability refers 
to the ability to be able to decode from the same 
stream to better and better qualities as more bits 
are decoded) in order to support a wide range of bit 
rates corresponding to available user bandwidths. 
The state-of-the-art image coder, 2D-SPIHT [16] 
produces a fully progressive code which means that 
if the transmission is stopped at any point, a lower 
bit image can still be decompressed and reconstructed. 
A longer prefix of the bit stream simply yields a 
higher reconstructed image quality when decoded. 
The reconstructed quality for such hierarchical bit- 
streams depends on which portions of the bitstream 
are received. In particular, if the first chunk of the 
stream is lost, then the remainder is essentially use- 
less. There is thus the need for an efficient mech- 
anism that converts a scalable, prioritized bitstream 
into a non-prioritized one so that is is better “matched” 
to the existing infrastructure. 

Multiple description (MD) source encoding has 
recently emerged as an attractive framework for ro- 
bust transmission over channels with transient shut- 
down characteristics. Examples of such channels 
are found in best-effort heterogeneous packet net- 
works such as the Internet, where congestion and 
routing delays can lead to “holes” in the transmit- 
ted data stream corresponding to lost packets, or 
wireless systems during a deep fade. For exam- 
ple, In the current deployment of packet networks, 
switches in the network are typically oblivious to 
the structure or content of the packets that they 
process and swatch. Packets that are more impor- 
tant to the application are treated identically to 
those that are not so important. In the following, 
w r e describe some recent w f ork that illustrates some 
key ideas for applying source-channel coding on an 
end-to-end basis. 

The basic idea in MD coding is to generate mul- 
tiple independent descriptions of the source such 
that each description independently describes the 
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source with a certain fidelity, and when more than 
one description is available, they can be synergis- 
tically combined to enhance the quality. The gen- 
eralized (n-channel) MD system can be informally 
characterized as a system that guarantees a hier- 
archy of delivered signal quality levels depending 
on how many descriptions k ^ n are correctly re- 
ceived as opposed to a traditional multiresolution 
(MR) system for which the quality delivered de- 
pends which k out of n layers are received. 

We consider a mechanism to transform a scal- 
able source bit stream into a robust MD packet 
stream by encoding source “layers” of decreasing 
importance with progressively weaker forward er- 
ror correction (FEC) [17] channel codes. Given the 
channel state and the transmission rate constraint, 
the question of how much protection to assign to 
each layer so as to maximize the expected quality at 
the receiver needs to be answered and we propose a 
nearly optimal solution to this problem. The major 
contribution of this approach is the rate-distortion 
“optimal” allocation between source and channel 
rates in the context of a best-effort packet network 
scenario. In the following, we describe the mechan- 
ics of the packetization strategy that converts the 
prioritized MR bitstream into an unprioritized MD 
stream using efficient erasure channel codes. 

The quality profile reflects the target quality (or 
equivalently distortion d: lower distortion implies 
higher quality and vice versa) when any k out of n 
descriptions are received. We will use the notation 
d(k)) to describe the quality profile where the ith 
entry in d(k) represents the target quality when i 
descriptions are received. 

Given n and d(fc), and a progressive bitstream, 
the stream is marked at n different positions (see 
Figure 8) which correspond to the attainment of the 
distortion levels d(k) 7 and is thus partitioned into 
n sections or resolution layers. The ith layer should 
be decodable when i descriptions get through i.e., 
the number of erasures does not exceed n — i. This 
is attainable with the use of the family of Reed- 
Solomon erasure-correction block codes which are 
characterized by the “optimal” code parameters (n, i , 
n — i + 1) which can correct any (n—i) erasures out 
of n descriptions. We split the i th quality layer into 

l An (n,k,d) block code is defined by a length n code 
with k user symbols and a minimum distance of d , i.e it can 
correct (rf— 1) erasures. Reed Solomon block codes have the 
property of maximum distance (d — n — k + 1) i.e the whole 
data can be recovered from any k out of n symbols. 



Figure 7: BER versus E^/Nq in a near- far setting 
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Figure 8: Progressive bitstream from the source coder 
partitioned into m layers or quality levels. 
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Figure 9: Forward Error Correction codes based MD 
codes for a quality profile of order n, i.e. having n 
quality levels or “sections.” 
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i equal parts, and apply the (n,i,n — i + 1) Reed- 
Solomon code to get the “contribution” from the ith 
level or section to each of the n descriptions. The 
“contributions” from each of the n quality levels or 
sections are concatenated to form the n descriptions 
(see Figure 9). Thus, every description contains all 
n layers, and all n descriptions are “equal” in infor- 
mation content as intended. 
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